Epigenetics & Chromatin

Post-translational modifications (PTMs) of histones are important epigenetic determinants and specific core histones PTMs correlate with functional chromatin states. However, despite linker histone H1s are heavily post-translationally modified, little is known about the genomic distribution of H1s PTMs and their association with epigenetic chromatin states. Here, we address this question in Drosophila that encodes a single somatic linker histone, dH1. We previously reported that dH1 is dimethylated at K27 (dH1K27me2). Here, we show that dH1K27me2 is a major PTM of Drosophila heterochromatin. At mitosis, dH1K27me2 accumulates at pericentromeric heterochromatin, while, in interphase cells, it is also detected at intercalary heterochromatin. ChIPseq experiments show that dH1K27me2 enriched regions cluster at both the assembled and unassembled heterochromatin regions of all four Drosophila chromosomes. More than 98% of the dH1K27me2 enriched regions map to heterochromatic repetitive DNA elements, including transposable elements, simple DNA repeats and satellite DNAs. We also show that dH1K27me2 is independent of H3K9 methylation, as it is equally detected in flies carrying a H3K9R mutation. Moreover, dH1K27me2 is not affected by depletion of Su(var)3-9, HP1a and Su(var)4-20. Altogether these results suggest that dH1K27me2 is a novel epigenetic mark of Drosophila heterochromatin that acts upstream of the major Su(var)3-9/HP1a pathway of heterochromatin formation.

RNA polymerase II (RNAPII) pausing at gene promoters is a rate-limiting step in transcription regulation. Previous studies have elucidated the coordinated actions of pausing and releasing factors that collectively modulate RNAPII pausing. In general, the involvement of chromatin remodellers in RNAPII pausing has not been well documented. Whilst LSD1 is well-known for its role in decommissioning enhancers during ESC differentiation, its role at the promoters of genes remains poorly understood despite their widespread presence at these sites. Here, we report that LSD1 is associated with RNAPII pausing at the promoter-proximal region of genes in mouse embryonic stem cells (ESCs). We demonstrate that the knockdown of LSD1 preferentially affects genes with higher RNAPII pausing than those with lower pausing and, importantly, show that the co-localization of LSD1 and MYC, a factor known to regulate pause-release, is associated with the enrichment of other RNAPII pausing factors compared to their independent counterparts. Moreover, we found that genes co-occupied by LSD1 and MYC are significantly enriched for housekeeping genes that are involved in metabolic processes and globally depleted of transcription factors compared to those bound only by LSD1. These findings reveal a pleiotropic role of LSD1 in regulating housekeeping program besides its previously known role in regulating cell identity programs. Our integrative analysis presents evidence for a previously unanticipated role of LSD1 in RNAPII pausing through its association with pause release factors in modulating cell-type specific and cell-type invariant genes.

BackgroundKDM5 family proteins are multi-domain regulators of transcription that when dysregulated contribute to cancer and intellectual disability. KDM5 proteins can regulate transcription through their histone demethylase activity in addition to demethylase-independent gene regulatory functions that remain less characterized. To expand our understanding of the mechanisms that contribute to KDM5-mediated transcription regulation, we used TurboID proximity labeling to identify KDM5-interacting proteins. ResultsUsing Drosophila melanogaster, we enriched for biotinylated proteins from KDM5-TurboID-expressing adult heads using a newly generated control for DNA-adjacent background in the form of dCas9:TurboID. Mass spectrometry analyses of biotinylated proteins identified both known and novel candidate KDM5 interactors, including members of the SWI/SNF and NURF chromatin remodeling complexes, the NSL complex, Mediator, and several insulator proteins. ConclusionsCombined, our data shed new light on potential demethylase-independent activities of KDM5. In the context of KDM5 dysregulation, these interactions may play key roles in the alteration of evolutionarily conserved transcriptional programs implicated in human disorders.

Core histone proteins H2A, H2B, H3, and H4 are encoded by a large family of genes distributed across the human genome. Canonical core histones contribute the majority of proteins to bulk chromatin packaging, and are encoded in 4 clusters by 65 coding genes comprising 17 for H2A, 18 for H2B, 15 for H3, and 15 for H4, along with at least 17 total pseudogenes. The canonical core histone genes display coding variation that gives rise to 11 H2A, 15 H2B, 4 H3, and 2 H4 unique protein isoforms. Although histone proteins are highly conserved overall, these isoforms represent a surprising and seldom recognised variation with amino acid identity as low as 77% between canonical histone proteins of the same type. The gene sequence and protein isoform diversity also exceeds commonly used subtype designations such as H2A.1 and H3.1, and exists in parallel with the well-known specialisation of variant histone proteins. RNA sequencing of histone transcripts shows evidence for differential expression of histone genes but the functional significance of this variation has not yet been investigated. To assist understanding of the implications of histone gene and protein diversity we have catalogued the entire human canonical core histone gene and protein complement. In order to organise this information in a robust, accessible, and accurate form, we applied software build automation tools to dynamically generate the canonical core histone repertoire based on current genome annotations and then to organise the information into a manuscript format. Automatically generated values are shown with a light grey background. Alongside recognition of the encoded protein diversity, this has led to multiple corrections to human histone annotations, reflecting the flux of the human genome as it is updated and enriched in reference databases. This dynamic manuscript approach is inspired by the aims of reproducible research and can be readily adapted to other gene families.

Acetylation of lysine residues in the tail domain of histone H3 is well characterized, but lysine residues in the histone globular domain are also acetylated. Histone modifications in globular domain have regulatory potential because of their impact on nucleosome stability but remain poorly characterized. In this study we report the genome-wide distribution of acetylated H3 lysine 115 (H3K115ac), a residue on the lateral surface at the nucleosome dyad, using chromatin immunoprecipitation. In mouse embryonic stem cells, we find that detectable H3K115ac is enriched at the transcription start site of active CpG island promoters, but also at polycomb repressed promoters prior to their subsequent activation during differentiation. By contrast, at enhancers H3K115ac enrichment is dynamic, changing in line with gene activation and chromatin accessibility during differentiation. Most strikingly, we show that H3K115ac is detected as enriched on "fragile" nucleosomes within nucleosome depleted regions at promoters, and active enhancers where it coincides with transcription factor binding, and at CTCF bound sites. These unique features suggest that H3K115ac correlates with, and could contribute, to nucleosome destabilization and that it might be a valuable marker for identifying functionally important regulatory elements in mammalian genomes.

BackgroundThe histone variant H3.3 is enriched at active regulatory elements such as promoters and enhancers in mammalian genomes. These regions are highly accessible, creating an environment that is permissive to transcription factor binding and the recruitment of transcriptional coactivators that establish a unique chromatin post-translational landscape. How H3.3 contributes to the establishment and function of chromatin states at these regions is poorly understood. ResultsWe performed genomic analyses of features associated with active promoter chromatin in mouse embryonic stem cells (ESCs) and found evidence of subtle yet widespread promoter dysregulation in the absence of H3.3. Loss of H3.3 deposition at promoters reduces chromatin accessibility and transcription factor (TF) footprinting for nearly all TFs expressed in ESCs. H3.3 deletion leads to reduced promoter enrichment of the transcriptional coactivator and histone acetyltransferase, p300. Subsequently, histone H3 acetylation at lysine 27 (H3K27ac) is reduced at promoters in the absence of H3.3, along with reduced enrichment of the bromodomain-containing protein BRD4, an acetyl lysine reader. Despite the observed chromatin dysregulation, H3.3 KO ESCs maintain transcription from ESC-specific genes. However, upon undirected differentiation, H3.3 KO cells retain footprinting of ESC-specific TFs motifs and fail to generate footprints of lineage-specific TF motifs, in line with their diminished capacity to differentiate. ConclusionsH3.3 facilitates DNA accessibility, TF binding, and histone post-translational modification at active promoters. While H3.3 is not required for maintaining transcription in ESCs, it is required for TF binding at new promoters during differentiation.

Histone H1 is involved in the regulation of chromatin structure and gene expression. Up to seven H1 subtypes or variants are expressed in human somatic cells. The H1 complement, defined as the subtype composition and proportions in a given cell, is highly variable depending on the cell type, cell cycle, developmental stage, and several diseases such as cancer. It is the result of the combined action of the different regulatory processes. Epitranscriptome modifications have emerged as a new regulatory mechanism able to control all aspects of mRNA metabolism. In this work, we have examined the role of the most prevalent mRNA modification, N6-methyladenosine RNA methylation (m6A), in the regulation of H1 subtypes. MeRIP-seq showed that H1.0 and H1.4 are enriched in m6A, while H1.2 has medium levels of this mark. We found that m6A inhibition altered transcript and protein levels, de novo transcription, and ribosome occupancy at the translation start site of specific H1 subtypes. Pull-down experiments using biotinylated-specific probes followed by mass spectrometry or RNA-immunoprecipitation coupled with RT-qPCR, allowed the identification of IGF2BP1, hnRNPD, and YTHDF2 as m6A readers of H1 subtypes. Integration of the functional studies, with m6A inhibition and partial depletion of the m6A readers led us to propose the first model of the differential regulation of H1 subtypes by m6A. In this model, m6A promotes the degradation of H1.0 mRNA mediated by YTHDF2, the stabilization of H1.2 mRNA by IGF2BP1 binding, and the H1.4 mRNA translation mediated by hnRNPD. These findings suggest that m6A is involved in the subtype-specific regulation of H1 subtypes, adding another layer to their complex regulation and contributing to the variability of the H1 complement in cancer.

We explored the role of NAT10 acetyltransferase in the DNA damage response, focusing on its impact on 3D-genome architecture and DNA repair proteins. Compared to NAT10 wild-type (wt), NAT10 deficiency reduced XPC, DDB2, and p53 protein levels. In TP53 double-null (dn) cells, the NAT10 protein was undetectable, and DDB2 was significantly down-regulated. Although NAT10 depletion caused DDB2 down-regulation, it did not affect the DNA repair functions of the DDB2 protein. To this fact, protein interaction analysis revealed that UVC exposure weakens the DDB2-p53 interaction while strengthening the bond between NAT10 and DDB2. Also, AlphaFold 3 prediction tools showed a more potent interaction between DDB2 and p53 than DDB2 and NAT10 proteins implying that NAT10 rather regulates the DDB2-p53 protein complex. These proteomic NAT10-dependent changes coincided with alterations in chromatin interactions, particularly in acrocentric chromosomes, studied by the Hi-C technique. However, 3D-genome rearrangement, caused by NAT10 deficiency and UVC irradiation, did not significantly impact post-translational histone modifications. Overall, NAT10 depletion alters the pool of key DNA repair proteins and induces substantial 3D-genome reorganization. Graphical abstractThe effect of the NAT10 acetyltransferase depletion and UVC irradiation on 3D-genome nuclear architecture, histone signature, and the pool of selected DNA repair-related proteins. The figure was made using some icons adapted from BioRender software.

BackgroundThe number and escape levels of genes that escape X chromosome inactivation (XCI) in female somatic cells vary among tissues and cell types, potentially contributing to specific sex differences. Here we investigate the role of CTCF, a master chromatin conformation regulator, in regulating escape from XCI. CTCF binding profiles and epigenetic features were systematically examined at constitutive and facultative escape genes using mouse allelic systems to distinguish the inactive X (Xi) and active X (Xa) chromosomes. ResultsWe found that escape genes are located inside domains flanked by convergent arrays of CTCF binding sites, consistent with the formation of loops. In addition, strong and divergent CTCF binding sites often located at the boundaries between escape genes and adjacent neighbors subject to XCI would help insulate domains. Facultative escapees show clear differences in CTCF binding dependent on their XCI status in specific cell types/tissues. Concordantly, deletion but not inversion of a CTCF binding site at the boundary between the facultative escape gene Car5b and its silent neighbor Siah1b resulted in loss of Car5b escape. Reduced CTCF binding and enrichment of a repressive mark over Car5b in cells with a boundary deletion indicated loss of looping and insulation. In mutant lines in which either the Xi-specific compact structure or its H3K27me3 enrichment was disrupted, escape genes showed an increase in gene expression and associated active marks, supporting the roles of the 3D Xi structure and heterochromatic marks in constraining levels of escape. ConclusionOur findings indicate that escape from XCI is modulated both by looping and insulation of chromatin via convergent arrays of CTCF binding sites and by compaction and epigenetic features of the surrounding heterochromatin.

Pericentromeric heterochromatin (PCH) is delineated by the enrichment of repressive epigenetic modifications, specifically trimethylated histone H3 at lysine 9 (H3K9me3) and DNA methylation (5-methylcytosine), which establish and maintain a condensed, transcriptionally silenced chromatin state. Depletion of either H3K9me3 or DNA methylation in mouse embryonic stem cells (mESCs) induces a permissive chromatin configuration that permits de novo recruitment and deposition of normally excluded Polycomb Repressive Complexes 1 and 2 (PRC1 and PRC2), characterized by H2AK119ub1 and H3K27me3 modifications, respectively, at PCH. Here, we demonstrate that H2AK119ub1 and H3K27me3 are independently recruited to hypomethylated PCH using a doxycycline-inducible mESC model allowing modulation of Dnmt1 expression levels and catalytic activity. We further investigate the roles of proposed mediators of PRC1/2 targeting, including SCML2, BEND3, KDM2b, and TET enzymes, in this context, our findings indicate that neither PRC1 nor PRC2 recruitment at hypomethylated PCH depends on these factors. Additionally, our data suggest that the permissive chromatin environment resulting from DNA hypomethylation is the principal facilitator of Polycomb complex spreading, offering novel insights into the mechanisms governing epigenetic modifier dynamics and interactions during periods of DNA methylation reprogramming.

Rubinstein-Taybi syndrome (RSTS) is an autosomal dominant disorder with specific clinical signs and neurodevelopmental impairment. The two known proteins altered in the majority of RSTS patients are the histone acetylation regulators CBP and p300. For assessing possible ameliorative effects of exogenous and endogenous HDAC inhibitors (HDACi), we exploited in vivo and in vitro RSTS models. First, HDACi effects were tested on Drosophila melanogaster, showing molecular rescue. In the same model, we observed a shift in gut microbiota composition. We then studied HDACi effects in RSTS cell lines compared to healthy donor cells. We observed patients-specific molecular rescue of acetylation defects at subtoxic concentrations. Finally, we assessed commensal gut microbiota composition in a cohort of RSTS patients compared to healthy siblings. Intriguingly, we observed a significant depletion in butyrate-producing bacteria in RSTS patients. In conclusion, this study reports the possibility of modulating acetylation equilibrium by HDACi treatments and the importance of microbiota composition in a chromatinopathy.

The six PCGF proteins (PCGF1-6) define the biochemical identity of Polycomb Repressor Complex 1 (PRC1) subcomplexes. While structural and functional studies of PRC1 subcomplexes have revealed specialized roles in distinct aspects of epigenetic regulation, our understanding of variation in protein interaction networks between the PCGF subunits is incomplete. We carried out an affinity purification mass spectrometry (AP-MS) screen of subunits PCGF1 (NSPC1), PCGF2 (MEL18), and PCGF4 (BMI1), using an immunoprecipitation approach that replicated endogenous cellular conditions in a cell line capable of differentiation programs. Over 200 interactions were found, including 83 that had not been described previously. Bioinformatic analysis found that these interacting proteins covered a range of functional pathways, often focused on cell biology and chromatin regulation. We found evidence of mutual regulation (at mRNA and protein level) between distinct PCGF subunits. Furthermore, we confirmed that disruption of each subunit using shRNA results in reduced proliferation ability. Overall, our work adds to understanding of the role of PCGF proteins within the wider cellular network.

The epigenetic modulator, DNMT3L, has been shown to be involved in nuclear reprogramming. Previous work from our laboratory had shown the accumulation and inheritance of epimutations across multiple generations when DNMT3L was ectopically expressed in transgenic DNMT3L Drosophila. Here, we show interaction of DNMT3L with Piwi, a member of a family of proteins known to perform their function through their association with piRNAs. Importantly, DNMT3L expression caused a significant alteration in the piRNA profile across multiple generations in transgenic Drosophila. As piRNAs are known to be passed on from one generation to another, we believe that the DNMT3L through its interaction with Piwi act to modify the inherited pool of piRNAs, which in turn allows accumulation of aberrant epigenetic modifications in subsequent generation. Furthermore, we show the interaction of DNMT3L with Histone H1, a non-core histone involved in higher order chromatin organisation. In light of these observations, it is proposed that in addition to its role in modulating core histone modifications, DNMT3L allows for inheritance of non-genetic information through its collaboration with Piwi, piRNAs and histone H1.

DNA methylation is extensively reprogrammed during early stage of mammalian development and is essential for normal embryogenesis. It is well established that mouse embryos acquire genome-wide DNA methylation during implantation, referred to as de novo DNA methylation, from globally hypomethylated blastocysts. However, the fact that the main de novo DNA methyltransferase 3B (DNMT3B) is initially expressed as early as the 8-cell stage, contradicts the current knowledge about timing of initiation of de novo DNA methylation. Here, we reported that a previously overlooked minor wave of de novo DNA methylation initially occurs during the transition from the 8-cell to blastocyst stage, before the well-known large-scale de novo DNA methylation during implantation. Functional analyses indicated that minor de novo DNA methylation regulates proliferation, lineage differentiation and metabolic homeostasis of preimplantation embryos, and is critical for embryonic developmental potential and pregnancy outcomes. Furthermore, bioinformatic and functional analyses indicated that minor de novo DNA methylation preferentially occurs on the X chromosome and co-regulates imprinted X-chromosome inactivation via the interaction between DNMT3B and polycomb repressive complexes 2 core components during blastocyst formation. Thus, our study updates the current knowledge of embryonic de novo DNA methylation, thereby providing a novel insight of early embryonic epigenetic reprogramming. Summary statementA minor wave of de novo DNA methylation has been initiated prior to blastocyst formation, but not during the implantation period, and co-regulates imprinted X-chromosome inactivation.

Polycomb Repressive Complex 2 (PRC2), an important histone modifier and epigenetic repressor, has been known to interact with RNA for almost two decades. In our previous publication (Long, Hwang et al. 2020), we presented data supporting the functional importance of RNA interaction in maintaining PRC2 occupancy on chromatin, using comprehensive approaches including an RNA-binding mutant of PRC2 and an rChIP-seq assay. Recently, concerns have been expressed regarding whether the RNA-binding mutant has impaired histone methyltransferase activity and whether the rChIP-seq assay can potentially generate artifacts. Here we provide new data that support a number of our original findings. First, we found the RNA-binding mutant to be fully capable of maintaining H3K27me3 levels in human induced pluripotent stem cells. The mutant had reduced methyltransferase activity in vitro, but only on some substrates at early time points. Second, we found that our rChIP-seq method gave consistent data across antibodies and cell lines. Third, we further optimized rChIP-seq by using lower concentrations of RNase A and incorporating a catalytically inactive mutant RNase A as a control, as well as using an alternative RNase (RNase T1). The EZH2 rChIP-seq results using the optimized protocols supported our original finding that RNA interaction contributes to the chromatin occupancy of PRC2.

The Illumina Infinium MethylationEPIC v2.0 BeadChip (EPICv2 array) is a microarray for assessment of the human epigenome. Sites on the EPICv2 array are annotated with an open-source file provided by Illumina, the EPICv2 manifest. Of the 923,452 unique genomic sites targeted by the EPICv2 array, the Illumina manifest identifies just 214,808 as mapping to a gene, excluding many sites located within a gene body. Based on the genomic coordinates of probes, we have mapped each site assayed on the Illumina EPICv2 array using publicly available data, comprehensively annotating affiliated genes and regulatory elements. We have found that a total of 700,392 EPICv2 array sites are located within a gene body (exon, intron, or UTR) according to the GENCODE Human release 47 (GENCODEv47) database. 509,940 of these sites were not annotated as being within a gene in the Illumina EPICv2 manifest, primarily because the Illumina manifest does not annotate introns - 498,407 of the excluded sites, or 97.74%, are located within the intron of at least one transcript. The Illumina EPICv2 manifest annotates 358,539 sites as being within 1500bp of a transcription start site (TSS). Using a distance-based approach, we have labelled 267,183 sites as being within promoter distance of a gene (<1500bp upstream or <500bp downstream of the TSS), and 140,123 sites as being within enhancer distance (1501-5000bp upstream of the TSS, excluding sites located within a gene body). We re-annotated the EPICv2 manifest using GENCODEv47 data to label intragenic features, and a distance-based approach to label the regulatory genome. We also include a column indicating whether a site is located in any promoter or enhancer, according to the GeneHancer database. The re-annotated manifest additionally labels which sites are required for the Horvath DNA Methylation Age Calculator and MethylDetectR epigenetic clocks, to facilitate data preparation for these tools. In conclusion, we have re-annotated the EPICv2 manifest, allowing more complete assessment of EPICv2 sites associated with gene bodies and regulatory regions during the interpretation of epigenetic studies. The re-annotated manifest is publicly available - see the Data Availability section of this article.

DNA methylation is one of the most important epigenetic processes that regulates gene expression. While the human genome is predominately methylated, CpG-rich regions known as CpG islands (CGIs) are unmethylated and are sites of transcriptional control. While the importance of CGI and DNA methylation at gene promoters is well understood, the significance of CGIs at 3 untranslated regions (3UTRs) remains largely unexplored. In this study, we characterised CGIs located within 3UTRs and investigated their role in gene regulation. Here, we show that around 3% of CGIs (909 CGIs) are exclusively located at 3UTR and are not associated with any nearby promoters. Importantly, 3UTR CGIs are highly conserved, with a conservation score even higher than promoter CGI implying their importance. In contrast to promoter CGIs which are predominantly unmethylated, 3UTR CGIs are often methylated. Genes with 3UTR CGIs are associated with several different cancers and cancer-related signalling pathways. Moreover, 3UTR CGIs are differentially methylated in cancers, with changes in methylation associated with changes in gene expression. Together this data suggests the importance of 3UTR CGIs and their methylation in gene expression regulation in cancer.

The SWI/SNF chromatin remodeling complex consists of more than 10 component proteins that form a large protein complex of > 1 MDa. The catalytic proteins Smarca4 or Smarca2 work in concert with the component proteins to form a chromatin platform suitable for transcriptional regulation. However, the mechanism by which each component protein works synergistically with the catalytic proteins remains largely unknown. Here, we report on the function of Smarce1, a component of the SWI/SNF complex, through the phenotypic analysis of homozygous mutant embryonic stem (ES) cells. Disruption of Smarce1 induced the dissociation of other complex components from the SWI/SNF complex. Histone binding to DNA was loosened in homozygous mutant ES cells, indicating that disruption of Smarce1 decreased nucleosome stability. Sucrose gradient sedimentation analysis suggested an ectopic genomic distribution of the SWI/SNF complex, accounting for the misregulation of chromatin conformations. Unstable nucleosomes remained during ES cell differentiation, impairing the heterochromatin formation that is characteristic of the differentiation process. These results suggest that Smarce1 guides the SWI/SNF complex to the appropriate genomic regions to generate chromatin structures adequate for transcriptional regulation.

RBBP4 is a core subunit of polycomb repressive complex 2 (PRC2) and HDAC1/2-containing complexes, which are responsible for histone H3 lysine 27 (H3K27) methylation and deacetylation respectively. However, the mechanisms by which RBBP4 modulates the functions of these complexes remain largely unknown. We generated viable mouse embryonic stem cell lines with RBBP4 mutations that disturbed methylation and acetylation of H3K27 on target chromatin and found that RBBP4 is required for PRC2 assembly and H3K27me3 establishment on target chromatin. Moreover, in the absence of EED and SUZ12, RBBP4 maintained chromatin binding on PRC2 loci, suggesting that the pre-existence of RBBP4 on nucleosomes serves to recruit PRC2 to restore H3K27me3 on newly synthesized histones. As such, disruption of RBBP4 function led to dramatic changes in transcriptional profiles. In spite of the PRC2 association, we found that transcriptional changes were more closely tied to the deregulation of H3K27ac rather than H3K27me3 where increased levels of H3K27ac were found on numerous cis-regulatory elements, especially putative enhancers. These data suggest that RBBP4 controls acetylation levels by adjusting the activity of HDAC complexes. As histone methylation and acetylation have been implicated in cancer and neural disease, RBBP4 could serve as a potential target for disease treatment.

Transcriptional silencers are cis-regulatory elements that downregulate the expression of target genes. Although thousands of silencers have been identified experimentally, a predictive chromatin signature of silencers has not been found. H4K20me1 previously was reported to be highly enriched among human silencers, but our reanalysis of those data using an appropriate background revealed that the enrichment is only marginal. We generated H4K20me1 ChIP-seq profiles in Drosophila S2 cells, which similarly showed that H4K20me1 does not mark Drosophila silencers and instead is associated with active transcription. Silencers remain a poorly annotated, difficult to predict class of cis-regulatory elements whose specific chromatin features remain to be identified.